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ABSTRACT 


The  heaving  and  pitching  motions  of  ALVIN  and  her  support  catamaran  are 
predicted  theoretically  and  experimentally  determined  at  zero  speed  by  model 
tests.   Comparison  is  made  between  theory  predictions  and  experimental  results 
for  each  of  the  vessels  and  then  the  model  test  results  for  the  two  vessels, 
with  ALVIN  in  the  recovery  position  are  compared. 

The  theoretical  results  were  computed  by  a  computer  program  based  on  the 
Korvin-Kroukovsky  linear  theory  of  ship  motions  in  conjunction  with  Grim's 
added  mass  and  damping  coefficients.   The  experimental  work  with  1/20  scale 
models  was  performed  at  the  M.I»T.  Ship  Model  Towing  Tank» 

It  is  concluded  that  theory  predictions,  correlates  reasonably  well  with 
model  test  results  for  the  catamaran.   The  comparisons  for  ALVIN  were  invali- 
dated by  questionable  experimental  data  and  possible  inapplicability  of  theory, 
From  experimental  tests,  it  was  concluded,  that  the  motions  of  the  two  ships 
are  synchronous  at  wavelengths  about  equal  to  the  length  of  the  catamaran, 
that  the  catamaran  significantly  damps  the  heave  of  ALVIN,  and  that  the  ALVIN 
has  little  effect  on  the  catamaran,  in  the  recovery  position.   Recommendations 
for  further  research  in  this  important  part  of  a  "Deep  Submergence  Vehicle 
System"  are  presented • 


Thesis  Supervisor:  Martin  A.  Abkowitz,  Ph°D< 
Title:   Professor  of  Naval  Architecture 
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CHAPTER  I 


INTRODUCTION 


The  deep  submergence  vehicle  is  but  a  component  of  a  larger  complex,  known 
as  the  "Deep  Submergence  Vehicle  System."  Next  to  the  submarine,  the  surface 
support  unit  is  the  most  important  component  of  this  system.   The  support  unit 
provides  transportation,  maintenance,  launching  and  recovery,  surface  tracking, 
navigation,  and  communications  with  the  submarine. 

The  launch  and  recovery  of  the  submarine  in  the  open  ocean  presents  a 
challenging  problem  to  the  designer  of  such  a  system.   The  essence  of  the  prob- 
lem is  how  to  handle  the  heavy  and  unwieldy  load  of  the  submarine  safely  and 
without  damage  in  a  seaway. 

Wood  Hole  Oceanographic  Institute  designed  and  built  a  catamaran  support 
vessel  (Figure  I)  during  the  construction  of  ALVIN,  their  deep  submergence 
vehicle  (Figure  II) .   The  catamaran  was  equipped  with  a  cradle  that  could  be 
lowered  and  raised  between  the  hulls  for  launch  and  recovery  of  ALVIN.   This 
was  thought  to  be  the  best  solution  to  handling  of  ALVIN  (Dry  Weight-13.6  long 
tons)  in  a  seaway. 

During  the  design  of  the  catamaran  several  studies  were  conducted  on 
methods  to  absorb  the  shock  expected  on  initial  contact  between  the  submarine 
and  catamaran  cradle .   Because  initial  trials  on  launching  and  recovering  a 
dummy  ALVIN  in  a  seaway,  revealed  that  the  expected  shock  did  not  exist,  these 
studies  were  not  implemented 

The  catamaran  concept  has  proven  itself  as  an  extremely  successful  support 
vessel  in  the  open  ocean.   It  was  able  to  launch  and  recover  ALVIN  in  up  to 
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state  4  seas  and  to  be  a  stable  and  seaworthy  vessel  in  25-ft  seas 

Two  main  objectives  comprise  the  scope  of  this  investigation.   The  first 
is  an  evaluation  of  the  accuracy  of  linearized  theory  of  heaving  and  pitching 
motions  of  ALVIN  and  the  catamaran  by  comparing  theoretical  computations  to 
experimental  model  motion  data  under  similar  conditions.   The  second  was  to 
compare  the  ALVIN  model  motion  data  with  the  catamaran  model  motion  data  in  the 
recovery  position. 

The  conclusions  from  this  investigation  may  prove  useful  in  the  future 
design  of  a  "Deep  Submergence  Vehicle  System." 

FIGURE   I 


CATAMARAN 

DSRVT-l 


L0A-96FEET 
BEAM-48  FEET 
DISPLACEMENT-350T0NS 
DRAFT,  KEEL- 6  FEET 
DRAFT,  PROP- 8  FT.  6  IN. 


CATAMARAN  OUTBOARD  PROFILE 


(APPENDAGES  WERE  NOT  MODELED) 


(Fig.  16  in  "ALVIN,  6000  FT.  Submergence  Research  Vehicle") 

by  Mavor»  etc. 
SNAME,  1966 


-2- 


FIGURE   II 


O  I  2  3*  5  6  7  8  9  10         II         \Z         13         14         15        16        17         II 


If  It         17         It         19        20 


ALVIN  HULL  LINES 


APPENDAGES  WERE  NOT  MODELED,  BUT  PROPELLOR 
SHROUD  AND  EXTERNAL  BOUYANCY  BLOCKS  (NOT 
SHOWN)  WERE  MODELED. 


(Fig.  8  in  "ALVIN,  6000  FT.  Submergence  Research  Vehicle") 
by  J.N.  Mavor,  Froehlich,  Marquet,  Rainnie 
SNAME,  1966 
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CHAPTER  II 


PROCEDURE 


General 

This  study  was  conducted  by  two  methods: 

1)  the  theoretical  computations  of  heaving  and  pitching  motions 
of  ALVIN  and  of  the  catamaran., 

2)  the  experimental  investigation  of  heaving  and  pitching  motions 
of  ALVIN  and  catamaran  models  in  the  M.LT.  Ship  Model  Towing 
Tank. 

Presented  in  this  section  are  the  various  assumptions  and  procedures 
utilized  in  the  theoretical  computations  and  the  tests  conducted  with  the 
models. 

Theoretical 

The  theoretical  computations  were  based  upon  linear  strip  theory 

(3) 
developed  by  Korvin-Kroukovsky  "   in  conjunction  with  Grim  s  results  on  added 

(If.) 
mass  and  damping    .   This  theory  has  been  utilized  in  a  digital  computer 

program  developed  at  the  M.LTa,  Department  of  Naval  Architecture  and  Marine 

(7) 
Engineering  by  Haslum  and  Vassilopoulos   «   The  analytical  details  of  this 

linear  theory  as  utilized  in  the  computer  program  are  summarized  in  Appendix 

B  from  reference  (7).   Use  of  this  computer  program,  described  in  Appendix 

C,  was  made  for  the  theoretical  computations  => 

The  basic  input  information  required  for  ships  to  the  computer  program 

is  the; 

1)   number  of  stations,  length,  and  displacement 
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2)  maximum  beam,  sectional  area  coefficient,  and  draft  for  each 
station 

3)  radius  of  gyration,  and  longitudinal  location  of  center  of 
gravity 

The  assumptions  made  in  order  to  submit  ALVIN  and  catamaran  input  data  were; 

1)  catamaran  represented  as  a  single  hull  with  half  the  displace- 
ment of  the  catamaran.   Equivalent  to  assuming  no  hull  inter- 
ference o 

2)  the  maximum  beam  for  both  ships  occurred  at  the  waterline  for 
both  ships 

3)  ALVIN 's  shrouded  propellor  was  represented  as  a  solid  disk 

4)  submerged  stations  immediately  forward  of  the  shroud  were 
entered  as  if  a  vertically  solid  up  to  the  waterline  from  the 
maximum  beam. 

The  input  information  for  waves  to  the  computer  program  was: 

1)  wave  amplitude  -0d25  fto,  2.5  fto  full  scale 

2)  shortest  wavelength  -»75  fto,  15  ft»  full  scale 

3)  longest  wavelength  -  10  fto,  200  fto  full  scale 

4)  increment  steps  from  shortest  to  longest  wavelength  -0.25  ft., 
5  ft.  full  scale 

5)  speed  of  ship  -0=0 

The  computer  output  was,  for  each  wavelength  to  shiplength  ratio,  the 
frequency  of  encounter,  the  heave  amplitude  (fto),  the  heave  amplitude  non- 
dimensionalized  with  wave  amplitude,  the  pitch  (degrees),  the  pitch  (radians) 
non-dimensionalized  with  the  maximum  wave  slope,  2tt/X  h  .   In  addition, 
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although  not  used  in  this  investigation,  the  relative  bow  motion,  the  cosine 
component,  the  sine  component,  and  phase  angle  were  outputs. 

This  investigation  was  conducted  for  zero  forward  speed  of  ALVIN  and  the 
catamaran  corresponding  to  the  launch  and  recovery  conditions  of  the  system. 

Experimental 

The  experimental  investigation  was  conducted  in  the  M.I.T.  Ship  Model 
Towing  Tank  equipped  with  a  hydraulic  wave  generating  system. 

The  tank  is  rectangular  in  cross  section,  108  ft.  long,  8  ft.  7  in.  wide 
with  a  normal  water  depth  of  A  f to  An  electrically  driven  carriage  is  used 
for  a  point  of  model  attachment  and  has  an  instrument  patch  panel  for  connec- 
tion to  the  recording  equipment.   A  beach  at  the  opposite  end  from  the  wave 
generating  equipment  essentially  eliminates  wave  reflection. 

The  waves  are  generated  by  the  motion  of  a  paddle  at  one  end  of  the  tank, 
which  rotates  about  an  axis  at  the  tank  bottom.   The  paddle  fills  the  full 
width  and  depth  of  the  tank.   The  paddle  is  powered  by  a  high  pressure 
hydraulic  ram  and  the  motion  is  controlled  by  an  electronically  controlled 
servo  valve.   For  generating  a  train  of  regular  waves  of  a  given  length,  a 
sinusoidal  signal  from  a  previously  prepared  magnetic  tape  is  the  input  to 
the  wavemaker  control  system  where  the  desired  wave  amplitude  is  set  by  a 
potentiometer  reducing  the  5  volt  peak-to-peak  sinusoidal  signal.  A  more 
complete  description  given  by  Pearlman  may  be  found  in  Reference  (4). 

The  model  dimensions  were  a  compromise  governed  by  the  following  require- 
ments: 

1)   To  allow  the  testing  together  of  ALVIN  and  catamaran,  the  model 
scales  had  to  be  equal 
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2)  It  was  desirable  to  make  the  catamaran  model  less  than  5  ft.  for 
transportation  and  handling 

3)  The  ALVIN  model  should  be  as  large  as  possible  to  accommodate 
the  weight  of  instrumentation 

The  scale  chosen  was  1/20  the  full  scale  size. 

The  full  scale  and  model  particulars  are  listed  in  Table  I. 

TABLE  I 


Full  Scale  and  Model  Particulars 


Scale  =  1/20 

Item 

Full  Scale 

Model 

c 

Length,  L,ft. 

96.0 

4.8 

A 

Beam(overall),  f  t° 

48.0 

4.8 

T 

Beam(each  hull),  ft. 

14o0 

0.7 

A 

Draft,  T,ft. 

6.0 

0.3 

M 

Displacement,  tons 

352.5  (SW) 

0.0436  (FW) 

A 

LongL  C.Gc  aft. 

R 

of  F.P.,ft. 

48.0 

2.4 

A 

Longlo  radius  of 

N 

gyration, ft. 

0.91 

A 

Length,  L,  ft. 

22o0 

1.1 

L 

Beam,f to 

8.25 

0.4125 

V 

Draft,  T,ft. 

7.07 

0.354 

I 

N 

Displacement,  tons 
Longl.  C.Gc  aft. 

Of  FoPe,  ft. 

13  o  6  (SW) 
8.54 

0.001658  (FW) 
0.00215   (FW) 
due  to  approx. 
in  actual  model 

0.469 

Longl.  radius  of 
gyration,  ft. 

0.242 
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The  instrumentation  used  to  measure  wave  height,  pitching  motion  ampli- 
tude, and  heav ding,  motion  amplitude  is  now  described. 

Wave  Height 

During  the  author's  use  of  the  tow  tank,  the  sonic  wave  probe,  which  is 
the  usual  means  of  measuring  wave  height,  was  out  of  commission.   The  wave 
height  was  measured  by  a  two  wire  resistance  probe  that  was  extended  to  half 
depth  in  calm  water.   The  two  wires,  spaced  approximately  1  in.  apart  with 
the  water  completing  the  circuit,  is  one  half  of  a  resistance  bridge  with  the 
other  half  of  the  bridge  in  the  carrier  preamplifier .   The  varying  water 
height,  as  a  wave  passes  the  probe,  changes  the  resistance  of  the  bridge 
providing  an  electric  signal  proportional  to  the  change  in  water  level.   The 
electrical  signal  is  the  input  to  the  oscillograph  recorder,  thus  providing 
a  continuous  record  of  the  wave  system.   The  probe  was  linear  for  ±1.5  in.  of 
water  level  change  and,  with  a  scale  of  10  mm/ in,  the  oscillograph  recorder 
had  an  accuracy  of  less  than  ±0<>25  mm  or  within  2.5%. 

Pitching  Motion  Amplitude 

The  pitch  angle  record  was  obtained  using  a  pitch  bearing  mounted  in 
the  pitch  bearing  block  which  is  rigidly  attached  to  the  model.  A  rotary 
variable  differential  transformer  coil  is  attached  to  the  bearing  block.   The 
core  is  attached  to  and  rotates  with  the  bearings   Thus  the  angular  motion  of 
the  model  about  the  point  of  attachement  generates  a  signal  in  the  coil  which 
is  proportional  to  the  pitch  angle  and  is  another  input  to  the  oscillograph 
recorder.   This  provides  a  continuous  record  of  the  pitching  motion  as  a 
function  of  time. 


Heaving  Motion  Amplitude 

The  vertical  oscillations  of  the  model  were  measured  by  a  linear  variable 
differential  transformer.   The  core  was  within  an  aluminum  rod,  which  was  attached 
to  the  model.   The  rod  passes  through  a  ball  bearing  guide,  which  was  attached  to 
the  carriage,  and  contains  a  three  winding  coil 

A  heave  rod  mounted  securely  in  the  pitch  bearing  and  passing  through 
a  ball  bearing  guide  attached  to  the  carriage,  served  the  purpose  of  allow- 
ing vertical  motion  only,  while  the  pitch  bearing  permitted  pitch  angular 
rotation  only.   For  the  catamaran  model;  the  heave  rod  was  square  in  cross 
section  as  was  its  guide,  to  restrain  the  model  from  other  motions,  espe- 
cially yaw,   For  the  ALVIN  model;  the  tubular  heave  rod  had,  at  its  upper  end, 
a  vane  that  rode  in  a  teflon  slot  on  a  teflon  bearing,  to  restrain  the  model 
from  other  motions .   The  linearsyn  core  rod  is  fixed  to  the  heave  rod  with  a 
plexiglass  bracket  and  thus  follows  the  vertical  motions  of  the  model .   The 
linearsyn  differential  transformer  generates  a  signal  proportional  to  the 
vertical  motion  of  the  model  and,  as  another  input  to  the  oscillograph  record- 
er, provides  a  continuous  record  of  the  heaving  motion. 

Figure  III  is  a  picture  of  the  catamaran  model  showing  the  instrumenta- 
tion and  Figure  IV  is  a  picture  of  the  ALVIN  model  showing  the  instrumenta- 
tion.  Table  II  lists  the  instrumentation  weights.   Figure XXV  is  a  typical 
oscillograph  recording  tape  that  was  used  for  data  collection. 

The  models  were  ballasted  to  the  correct  draft  with  the  instrument 
weights.   Then  the  ballast  was  shifted  to  adjust  for  the  longitudinal  radius 
of  gyration.   The  procedure  used  to  determine  the  model  radius  of  gyration  is 
described  in  Appendix  D. 
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FIGURE  III 
CATAMARAN  MODEL  INSTRUMENTATION 


FIGURE  IV 
ALVIN  MODEL  INSTRUMENTATION 
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TABLE  II 


Instrumentation  Weights 


Model 

Item                Weight  (lbs) 

Catamaran 

Heave  Rod  (Solid) 

2.36 

Pitch  bearing  block,  bearing, 

.71 

and  mounting  bolts 

linearsyn  core  rod  and 

.385 

bracket 

Total 

3.455 

ALVIN 

Heave  Rod  (Hollow) 

.105 

Pitch  heaving  block,  bearing, 

.71 

and  mounting  bolts 

linearsyn  core  rod  and 

.385 

bracket 

Total 

1.19 

Test  Procedure 

To  determine  the  heave  and  pitch  response  of  ALVIN  and  the  catamaran, 
they  were  first  tested  separately  and  then  both  models  together  with  ALVIN 
in  the  recovery  position  between  the  catamaran  hulls.   The  runs  that  were 
conducted  are  listed  in  Table  III. 

All  runs  were  conducted  with  the  carriage  30  ft.  away  from  the  wavemaker 
paddle  to  ensure  regular  waves.   With  the  calibrated  wave  height  probe,  on 
the  carriage,  and  no  models  in  the  water,  the  wave  amplitude  settings  on  the 
potentiometer  were  determined  for  each  wavelength. 
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The  catamaran  model  was  tested  for  two  different  wave  amplitudes  at  each 
wavelength o   A  single  catamaran  hull  was  then  tested  in  the  same  manner. 
Because  of  the  bow  and  stern  symmetry  of  the  catamaran,  all  runs  were  for 
both  directly  ahead  and  astern  seas.   The  ALVIN  was  tested  in  directly  ahead 
and  directly  astern  seas. 

For   testing  together,   both  models  were  attached  to  the  carriage  with 
the  heave  rod  and  pitch  bearing.   This  allowed  both  models  to  move  freely  in 
vertical  motion  and  in  pitch  angular  rotation.   The  first  set  of  runs  were 
with  heave  and  pitch  instrumentation  on  ALVIN.   The  second  set  of  runs  were 
with  the  heave  instrumentation  on  the  catamaran  while  the  pitch  instrumenta- 
tion on  ALVIN.   This  method  of  attachment  is  shown  in  Figures  V  and  VI. 

TABLE  III 


Experimental  Runs 


Model        Wave  Amplitude  (ft.) 

Wavelength 

(ft.) 

Catamaran 

0.0625 

3.0, 

4.0,  4.5, 

6.0,  10.0 

0.125 

3.0, 

4.0,  4.5, 

6.0,  10.0 

Catamaran 

0.0625 

3.0, 

4.0,  4.5, 

6.0,  10.0 

(Single  hull) 

0.125 

ALVIN,  ahead 
seas 

0.0625 

1.0, 

2.0,  3.0, 

5.0,  10.0 

ALVIN,  astern 
seas 

0.0625 

1.0, 

2.0,  3.0, 

5.0,  10.0 

ALVIN,  Recovery 

Position,  ahead 

0.0625 

1.0, 

2.0,  3.0, 

5.0,  10.0 

seas 

ALVIN ,  Recovery 

Position,  ahead 

0.0625 

3.0, 

4.0,  4.5, 

6.0,  10.0 

seas,  pitching 

only 

Catamaran,  Recovery 
heaving  only 

0.0625 

3.0, 

4.0,  4.5, 

6.0,  10.0 
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FIGURE  V 
RECOVERY  POSITION  MODEL  INSTRUMENTATION 


FIGURE  VI 
RECOVERY  POSITION  MODEL  IN  WAVES 
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CHAPTER  III 
ANALYSIS 

The  theoretical  computed  results  for  the  catamaran  were  for  a  single  hull 
only,,   There  is  no  known  way  for  the  Korvin-Krovskovsky  linear  strip  theory  to 
take  into  account  two  hulls  as  in  a  catamaran  configuration.   Either  the  assump- 
tion must  be  made  that  no  interference  exists  between  the  two  hulls  or  a  correc- 
tion must  be  made. 

The  damping  of  heaving  and  pitching  motion  is  caused  by  the  oscillating 
body  generating  a  traveling  wave  system  in  which  energy  is  dissipated  (2). 
With  two  hulls  in  close  proximity  to  each  other  and  both  oscillating,  there 
must  be  interference  between  the  hulls.   Therefore  it  would  be  erroneous  to 
make  the  assumption  of  no  interference  and  accordingly  a  way  to  correct  the 
theoretical  results  was  sought. 

Because  there  was  no  theory  available,  empirical  methods  were  the  only 
means  available  and  the  procedure  used  to  correct  the  theoretical  results  is 
described  herein 

The  non-dimensional  experimental  heave  and  pitch  motion  amplitudes  for  the 
catamaran  and  a  single  hull  were  plotted  for  two  different  wave  heights  , 
(Figures  VII,  VIII,  IX,  and  X. )  From  the  data  obtained  for  the  two  different  wave 
heights,  at  each  wavelength,  a  mean  value  was  obtained  by  interpolation  for 
both  the  case  of  the  single  hull  and  the  catamaran.   These  points  are  plotted 
in  Figures  XI  and  XII .   The  ratio  of  the  catamaran  to  the  single  hull  was 
computed  for  the  different  wavelengths  and  applied  as  a  correction  to  the 
theoretical  non-dimensional  motion  amplitudes  for  the  single  hull.   (See 
Table  IV). 
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TABLE  IV 


Correction  for  Single  Hull  to  Catamaran 


Mean 

Mean 

Cat. 

Corrected 

X(ft) 

Single 

Cat. 

Single 

Theory 

Amplitudes 

(heave 

amplitude) 

3.0 

0 

,1200 

+.1200 

.2611 

.3811 

3,5 

o0700 

.0750 

1.072 

.2553 

.2735 

4.0 

.1034 

.0667 

0.645 

.1488 

.0960 

4.5 

o0200 

.0917 

4.58 

.0169 

.0775 

5.0 

0O8OO 

,1300 

1.625 

.1115 

.1815 

6.0 

.2083 

.2500 

1.2 

.3244 

.3895 

7.0 

.3000 

.3700 

1.232 

.4769 

.5880 

8.0 

.3800 

.5050 

1.33 

.5914 

.7860 

9.0 

»4425 

.6500 

1.47 

.6729 

.9880 

10.0 

,5000 

.8000 

1.6 

.7335 

1.172 

(pitch 

amplitude) 

3.0 

.0300 

.0501 

1.67 

.0499 

.0749 

3.5 

.0650 

.1600 

2.42 

.0536 

.1300 

4.0 

.1200 

.2245 

1.87 

.1645 

.3080 

4.5 

.2025 

.2755 

1.31 

.2689 

.3525 

5o0 

.2750 

.3625 

1.32 

.3608 

.4760 

6.0 

.375 

.4687 

1.25 

.5060 

.6325 

7.0 

,4525 

.5550 

1.228 

.6098 

.7480 

8.0 

.5250 

.6475 

1.232 

.6848 

.8430 

9.0 

.5950 

.7400 

1.242 

.7401 

.9200 

10.0 

.6590 

.8320 

1.262 

.7818 

.9850 
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CHAPTER  IV 


RESULTS 


The  results  are  presented  in  the  form  of  non-dimensional  motion  amplitudes 

versus  wavelength  and  wavelength  to  shiplength  ratio,  Figures  XIII-XXIVo   T1  've, 

z  ,  is  non-dimensionalized  with  the  wave  height,  h  .  and  pitch,  0  in  radians, 
o  o  o 

by  the  maximum  waveslope,  2tt/A  h  .   Full  lines  represent  theoretically  computed 
motion  amplitudes  and  experimentally  determined  motion  amplitudes  are  represented 
by  circled  points. 

Figures  XIII  and  XIV  compare  catamaran  theory  predictions  with  experimental 
results.   Figures  XV  and  XVI  show  ALVIN  comparisons  in  HeaVe  and  Figures  XVII  and 
XVIII  show  ALVIN  comparisons  in  Pitch.   Figure  XIX  is  a  comparison  of  ALVIN  with 
and  without  heave  instrumentation.   Figures  XX  and  XXI  are  comparisons  of  the 
two  vessels  in  the  recovery  position.   Figures  XXII  and  XXIII  are  comparisons  of 
ALVIN  alone  and  in  recovery  position.   Figure  XXIV  compares  the  catamaran  alone 
and  with  ALVIN  in  recovery  position  for  heave  motion. 
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CHAPTER  V 
DISCUSSION  OF  RESULTS 

For  the  catamaran,  after  the  theoretically  computed  heaving  motion 
amplitudes  had  been  corrected,  the  plotted  curve  correctly  predicted  the  trend 
of  experimental  results .  While  theory  predicted  somewhat  higher  at  shorter 
and  longer  wavelengths,  in  the  range  of  0.9  to  1.5  wavelength  to  ship  length 
ratio  there  was  reasonably  close  agreement  with  experimental  results    (Figure 
XIII), 

The  curve  of  corrected  theoretical  pitching  amplitudes  correctly  predicted 
the  trend  but  was  about  20%  higher  than  experimentally  determined  results . 
(Figure  XIV). 

Instrumenting  the  small  ALVIN  model  was  extremely  difficult  because  the 
weight  of  instrumentation  was  of  the  order  of  magnitude  of  weight  of  the  model. 
A  hollow  heave  rod  of  tubular  aluminum  was  substituted  for  the  solid  square 
steel  heave  rodo  However,  during  the  model  testing  for  motion  response,  the 
light  weight  of  the  instrumented  ALVIN  was  not  sufficient  to  overcome  the  small 
friction  in  the  guides  of  the  heave  rod  and  linearsyn  core  rod.   For  many  runs 
no  vertical  oscillation  was  recorded  at  all.   The  runs  where  heave  motion  was 
recorded  showed  intermediate  binding  so  that  all  experimental  heave  data  of 
ALVIN  was  considered  questionable . 

Because  of  the  undesired  restraint  on  vertical  motion,  the  model  responded 
to  the  wave  action  with  amplified  pitch  motion.  An  example  of  this  is  shown 
in  Figure  XIX,  which  compares  the  pitching  motion  with  and  without  the  linearsyn 
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rod.   This  proves  that  pitching  motion  was  artificially  increased  due  to  this 
unwanted  restraints   Another  piece  of  evidence  is  the  fact  that  at  5.0  ft,  wave- 
lengths and  10o0  ftc  wavelengths  pitch  exceeded  by  a  factor  of  2-5  the  maximum 
waveslope.   At  pitch  angles  about  equal  to  5  degrees  bow  up  the  center  of 
buoyancy  shift  was  great  enough  to  cause  an  instability  in  the  model  and  it 
would  flip  stern  down,,   For  this  reason,  the  maximum  wave  amplitude  that  could 
be  used  with  ALVIN  was  0.75  in,  crest  to  trough. 

Comparisons  with  theoretically  computed  heave  and  pitch  motions  amplitudes 
were  nevertheless  made,  albeit  somewhat  questionable.   Theoretical  results  were 
higher  for  both  directly  ahead  and  astern  seas.   Also  the  curve  of  theoretical 
results  of  heaving  amplitudes  predicted  a  resonance  at  the  2,25  ft.  wavelength 
that  was  not  verified  by  the  experimental  data,   (Figures  XV  and  XVI) .   For 
ALVIN  pitching  motions  the  experimentally  determined  results  crossed  the 
theoretical  curve  at  about  4.0  ft,  wavelength  in  both  directly  ahead  and  astern 
seas  where  the  restraint  in  vertical  motion  caused  greater  pitching  motion  of 
the  model,   (Figures  XVII  and  XVIII), 

The  recovery  position  motion  data  shows  that  for  heaving  amplitudes  the 
two  vessels  first  follow  the  same  trend  with  ALVIN  heaving  slightly  more  than 
the  catamaran  at  wavelengths  shorter  than  the  catamaran  length  then  about  the 
same  at  a  wavelength  90%  of  the  catamaran  length.   Then  the  curves  of  heaving 
amplitude  diverge  with  the  catamaran  heaving  amplitude  greater.   (Figure  XX). 
This  is  contrary  to  the  author's  observations  during  the  runs,  when  it  was 
directly  observed  at  the  4.5  ft.  wavelength  the  ships  seemed  to  move  in  perfect 
synchronism  and  exhibited  similar,  though  less  perfect,  synchronism  as  the  wave- 
length increased.   The  discrepancy  between  measured  results  and  observation  may 
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be  caused  by  binding  of  the  heave  instrumentation  on  the  ALVIN  model. 

For  the  recovery  position  pitching  motion,  the  rapid  divergence  of  ALVIN 
at  longer  wavelengths  is  felt  to  be,  again,  due  to  the  restraint  in  heave. 
(Figure  XXI) . 

ALVIN 's  effect  on  catamaran  heaving  motion  appears  to  be  slight  and 
highly  frequency  (wavelength)  dependent  <,      (Figure  XXII). 

ALVIN' s  heaving  motion  amplitude  appears  to  be  significantly  damped 
(Figure  XXIII). 

ALVIN 's  pitching  motion,  however,  is  amplified  at  the  shorter  wavelengths 
with  a  resonance  at  a  wavelength  62%  of  the  catamaran,  and  damped  slightly 
at  a  wavelength  90%  of  the  catamaran  length.   (Figure  XIV).   This  is,  again, 
inconclusive  because  of  the  ALVIN  model  instrumentation. 
Reasons  for  Discrepancies 

The  experimental  methods  and  the  current  status  of  theoretical  knowledge 
possess  many  sources  of  errors.   To  provide  a  partial  explanation  of  the 
differences  in  this  investigation  between  theory  and  experiment  and  the 
differences  between  the  models,  these  discrepancies  are  discussed. 

For  the  experimental  work,  the  following  sources  of  discrepancies  exist: 

a)  Wave  height  measurement  is  one  of  the  greatest  sources  of  errors 

(7) 
in  tow  tank  work   „   The  potentiometer  settings  for  wave- 
heights  were  determined  with  no  models  in  the  tank  to 
eliminate  the  model  effect  on  the  wave  height  measured.   How- 
ever, as  the  potentiometer  settings  were  quite  delicate,  this 
may  have  led  to  errors. 
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b)  Wall  effects,  or  reflection  from  the  walls  back  onto  the  model 

of  waves  generated  by  the  model  motion  at  zero  forward  speed 

(2  7) 

and  at  small  wavelengths  significantly  affected  the  results 

c)  Undesired  instrument  effects  on  the  model  may  have  affected  the 
results,  as  was  previously  discussed  for  the  ALVIN  models 

d)  Instrument  accuracy  was  about  ±2.5%. 

e)  Extraction  of  recorded  data  on  oscillograph  recorder  tapes  may 
have  introduced  errors.. 

f)  During  the  1.0  ft.  wavelength  generation,  the  hydraulic  ram 
pushing   the  paddle  was  observed  to  loosen  on  its  foundation 
and  a  modulated  wave  record  resulted.   This  may  have  affected 
other  wavelengths,  but  oscillograph  records  appeared  to  be 
regular o 

For  the  theoretical  computations,  the  following  assumptions  limited  the 
applicability  of  the  theory  to  ALVIN  and  the  catamaran: 

a)  linearized  theory  assumes  ship  section  lines  vertical  at  the 
waterline.  Neither  of  the  models  have  vertical  sides,  espe- 
cially in  the  case  of  ALVIN . 

b)  the  computation  of  damping  in  linear  theory  is  by  accounting 
for  the  energy  dissipated  in  gravity  waves  generated  when  the 
ship  oscillates „   Frictional,  separation,  and  circulation 

phenomena  are  assumed  to  be  negligible,  which  is  generally  valid 

(5) 
for  surface  ships   .   This  assumption  is  erroneous  for  ALVIN 

where  the  above  phenomena  are  a  considerable  influence. 
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c)  linear  theory  assumes  small  motions  and  that  coefficients  of 
motion  are  constant  throughout  the  whole  cycle  of  motion <, 
This  is  not  true  for  larger  motions.   The  restoring  force 
coefficients  are  correct  only  for  vertical  sides  at  the  water- 
line  o   For  the  ALVIN  model  with  a  decreasing  and  increasing 
waterplane  area  as  the  model  oscillates  up  and  down,  this  is 
certainly  not  true. 

d)  end  effects  or  three-dimensional  effects  are  responsible  for 
large  discrepancies  when  large  section  areas  are  present  at  the 
ends  of  the  model .   The  sectional  area  curve  must  go  "grace- 
fully" to  zero  at  the  bow  and  stern  in  order  for  the  force 
contributions  of  these  sections  to  zero.   This  does  not  affect 
the  heave  coefficients  seriously  but  the  ends  of  the  ship  have 

a  double  effect  on  pitching  and  therefore  significant  errors 

(5) 
may  result  in  pitching  motions   .   This  is  serious  inadequacy 

in  applicability  of  the  theory  to  ALVIN. 

e)  linear  strip  theory  has  a  basic  assumption  that  beam  and  draft 

are  small  compared  to  length  and  uses  a  basically  two-dimen- 

(5) 

sional  approach   .   This  is  another  serious  inadequacy  in  the 

applicability  of  the  theory  to  ALVIN,  which  has  a  length  to 
beam  ratio  of  2.6  and  a  length  to  draft  ratio  of  3.L 
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CHAPTER  VI 
CONCLUSIONS 

The  conclusions  from  the  preceding  chapters  are  summarized. 

1.  The  catamaran  corrected  theoretical  results,  in  spite  of  discrep- 
ancies noted  previously,  predicts  and  agrees  reasonably  well  with 
the  experimental  results. 

2.  The  Korvin-Krovkorsky  linear  theory,  with  Grim's  added  mass  and 
damping,  applicability  to  ALVIN  is  questionable  but  this  conclusion 
cannot  be  substantiated  because  of  inaccurate  experimental  data. 

3.  Because  of  the  inapplicability  of  theory  and  the  questionable 
experimental  data,  a  valid  comparison  could  not  be  made  between 
ALVIN  theoretical  predictions  and  experimental  results. 

4.  ALVIN  and  the  catamaran  respond  in  synchronous  motion  at  wavelengths 
about  equal  to  the  length  of  the  catamaran,  but  this  conclusion  is 
questionable  because  of  ALVIN  experimental  data. 

5.  ALVIN  heaving  motion  is  appreciably  damped  in  the  recovery  position 
between  the  catamaran  hulls,  but  this  conclusion  is  questionable 
because  of  ALVIN  experimental  data» 

6.  The  catamaran  heaving  motion  is  not  significantly  affected  by  ALVIN 
in  the  recovery  position . 
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CHAPTER  VII 


RECOMMENDATIONS 


Much  research  remains  to  be  done  in  the  design  of  surface  support  vessels 
for  lauching  and  recovery  of  deep  submersibles  in  the  open  ocean.   Therefore  it 
seems  important  that  further  theoretical  and  experimental  work  be  undertaken  in 
this  important  part  of  the  "Deep  Submergence  Vehicle  System." 

Recommendations  arising  from  this  investigation  are  that; 

1)  the  theory  be  modified  to  extend  its  applicability  to  deep  sub- 
mersibles in  surfaced  condition  similar  to  the  hull  form  of  ALVIN. 

(18) 
There  exists  in  the  literature  the  more  general  analysis  of  Porter 

/ION 

and  Bermejo     that  accounts  for  section  shape  in  more  detail. 

2)  the  comparisons  attempted  in  this  investigation  be  repeated  with 
special  attention  given  to  devising  better  heave  instrumentation 
for  ALVIN  or  larger  scale  models  be  used. 

3)  this  investigation  should  be  extended  by  experimentally  determining 
phase  angles  of  motion  and  comparing  them  with  the  theoretical 
computed  phase  angles  (that  were  not  used  in  this  investigation) . 

4)  an  experimental  investigation  be  made  using  different  support  vessel 
hull  forms  to  determine  the  optimum  hull  form  for  open  ocean. 
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APPENDIX  A 
NOMENCLATURE 

a  =  coefficient  of  equation  of  motion 

A  =  coefficient  of  equation  of  motion 

b  =  coefficient  of  equation  of  motion 

B  =  coefficient  of  equation  of  motion 

c  =  coefficient  of  equation  of  motion 

C  =  coefficient  of  equation  of  motion 

d  =  coefficient  of  equation  of  motion 

D  =  coefficient  of  equation  of  motion 

e  ■  coefficient  of  equation  of  motion 

E  =  coefficient  of  equation  of  motion 

F  =  total  vertical  (heaving)  force 

F  =  amplitude  of  time-varying  heaving  force 

F  =  complex  vertical  (heaving)  force 

g  =  gravitational  acceleration 

g  =  coefficient  of  equation  of  motion 

G  =  coefficient  of  equation  of  motion 

h  =  amplitude  of  sinusoidal  wave  (half-wave  height) 

K  =  maximum  wave  slope,  (2tt/X  h  ) 

K  =  longitudinal  radius  of  gyration 

o 

L  =  length  of  ship  or  model 

m  ■  mass  of  ship  or  model 

M  ■  total  pitching  moment 
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M      =  amplitude  of  time-varying  pitching  moment 

M      ■  complex  pitching  moment 

N(£)   ■  sectional  damping  coefficient 

t      -  time 

v      =  speed  of  ship  or  model 

z      =  heaving  motion  of  center  of  gravity  of  ship  or  model 

z      =  heaving  velocity  of  center  of  gravity  of  ship  or  model 

z°      =  heaving  acceleration  of  center  of  gravity  of  ship  or  model 
z      =  amplitude  of  heaving  motion 

z      ■  complex  heaving  motion 

Greek  Letters 

6  ■  theoretically  computed  heaving  phase  angle 

A  =  displacement  of  ship  or  model 

e  =  theoretically  computed  pitching  phase  angle 

6  =  pitching  displacement 

0  ■  pitching  velocity 

0  =  pitching  accelerations 

0  =  amplitude  of  pitching  motion 

0  =  complex  pitching  motion 

X  =     wavelength 

p  =  water  density 

O  =  phase  angle  of  heaving  force 

T  ■  phase  angle  of  pitching  moment 


w   =  frequency  of  encounter 
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APPENDIX  B 
ANALYTICAL  DETAILS  OF  THE  JLINEAR  THEORY  OF  SHIP  MOTIONS 

In  order  to  achieve  completeness  and  simultaneously  provide  easy  reference, 
the  fundamental  analytical  expressions  of  the  Korvin-Kroukovsky  method,  as 
utilized  in  the  computer  program,  are  summarized  in  this  Appendix.   Following 
the  nomenclature  and  definitions  adopted  in  (8) ,  the  coupled  set  of  linear 
differential  equations  describing  the  two-degree-of-freedom  ship  system,  takes 
the  form, 

az  +  bz  +  cz  +  d0  +  e0  +  g0  =  F  exp  (iw  t) 


AO  +  B0  +  C0  +  Dz  +  Ez  +  Gz  ■  M  exp  (iu)  t) 


(1) 


The  above  equations  result  from  equilibrium  considerations  of  the  hydro- 
dynamic  forces  and  moments  called  into  play  by  the  ship's  oscillations  in  the 
plane  of  symmetry,  when  meeting  head  or  astern  regular  waves.   For  this  reason 
the  analysis  ignores  steady,  continuously  acting  forces  due  to  buoyancy,  gravity 
and  suction  pressures.   Following  the  principles  of  classical  dynamics,  these 
forces  and  moments  are  obtained  by  applying  Newton's  Second  Law  of  Motion  to 
both  translatory  and  rotational  displacements  of  the  body's  center  of  gravity. 

The  wave  induced  excitation  force  and  moment  may  be  thought  of  as  being 
imposed  on  a  fully  restrained  ship  and  appear  on  the  RHS  of  the  set  (1). 
They  have  the  useful  property  that  they  are  functions  of  the  wave  elevation  and 
its  two  first  time  derivatives,  thereby  allowing  ease  in  algebraic  manipulation 
They  are  defined  in  complex  notation  as 


F  exp  (ico  t)  =  F  exp  (-ia)  exp  (ioo  t)  ■  F  exp  [i(oo  t-a)] 


M  exp  (ico  t)  =  Mq  exp  (-il)  exp  (ico  t)  =  M  exp  [i(w  t-i)] 
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(2) 


The  differential  exciting  force  acting  on  a  control  section  distant  £  from 
the  origin  of  the  moving  coordinate  system  (ship's  C.G.)»  has  been  given  in 
(8)  in  the  simplified  form, 


dF     dFl       r  dF2   .    , 

—  ■  - —  cos  a)  t  +  - —  sin  w  t 
dx     dx        e      dx        e 


=      [{<t>1  sin  -^    +  <f> 


— r —  cos  ~~y~s  exp   (-  — r'-Jcos  u  t 


_    r  27rhc 

^£     -  d>       S  sin  ^ 

'1  —     X  92  X       sin     X 


+     [{<$>.   cos  ^r5*     -  <J)0     — r-2,  sin  -^Hr3-}  exp    (-  -£rx)]sin  oo  t 


2& 
X 


(3) 


where, 


2        2 
4tt  he 

<|)1     =     hpg  B*  -     f~     (p  Sk2kA) 


and, 


<J>2     =     N(£)    -  V 


d(p  sk2k4) 

dC 


while   the  differential  exciting  moment  of   this   force  about   the  C.G.    is   given  by 

dF 

—  d£«   Integration  of  the  above  two  quantities  over  the  ship  length  results  in 

the  values  of  the  total  time-varying  exciting  force  and  moment,  which  are  con- 
sidered as  the  real  parts  of  the  expressions  (2).   Thus, 


F  =  F,  cos  0)  t  +  F_  sin  oj  t 
1      e     2      e 


M  =  M.  cos  co  t  +  M„  sin  0)  t 
1      e     2      e 


2     2  F 

F   +  F  2 

1     2  cos[oj  t-arctan  —  ]    and 

e        Fl 


M, 


2     2  2 

=  ,/M..   +  M_  cos[oj  t-arctan  — ] 
2       e        M.. 


F  cos(oj  t-a) 
o      e 


=  M  cos(w  t-x) 
o      e 


(A) 


-49- 


The  analysis  of  the  forces  and  moments  which  correspond  to  the  ship's  free 
oscillations  in  calm  water  yields  terms  which  appear  on  the  LHS  of  the  set  (1) 
and  are  proportional  to  the  instantaneous  heaving  and  pitching  displacement, 
velocities  and  accelerations.   All  twelve  coefficients  of  the  above  terms  are 
independent  of  the  speed  per  se  and  the  body's  space  orientation,  but  depend 
on  the  frequency  of  encounter,  with  the  exception  of  c  and  G.   The  final 
expressions  for  the  coefficients  of  the  equation  of  motion  used  in  the  machine 
computations  are  listed  below: 

a  =  m  +  p  /  (Sk2k4)  d£ 

b  =  /  N(C)  d5 

c  =  p  g  /  B*  d£ 

d  =  p  /  (Sk2k4)  E,   d£ 

e  =  /  N(£)  K   d£  -  2V  p  /  (Sk2k4)  d£  -  V  p  /  d  (Sk  k^/d^  E,   d£ 

g  =  P  g  /  B*  U^  -  V  /  N(£)  dE,  (5) 

A  =  J  +  p  /  (Sk2k4)  e}   di 

B  -  /  N(0  C2  dC  -  2V  p  /  (Sk2k4)  E,   d£  -  V  p  /  d  (Sk2k4)/d£  e}   d£ 

C  =   p  g  /  B*  e}   dE,   -  V  /  N(C)  5  d£  +  V2  p  /  d  (Sk2k4)/d£  E,   dC 

D  =  p  /  (Sk2k4)  5  d? 

E  =  /  N(£)  £  d£  -  V  p  /  d  (Sk2k4)/d£  C  d£ 
G  =   p  g  /  B*  £  d£ 

Assuming  now  that  sufficient  time  has  elapsed  for  any  transient  distur- 
bances to  be  damped  out,  we  seek  particular  solutions  of  the  nonhomogeneous  set 

(1),  which  correspond  to  the  steady-state  responses  of  the  system.   Since  the 
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forcing  functions  are  sinusoidal  and  the  system  is  linear  and  time-invariant, 
we  expect  that  any  response  will  also  be  sinusoidal  of  the  same  frequency  as 
the  excitation  and  with  generally  different  amplitude  and  phase.  We  therefore 
assume  solutions  of  the  form, 


(6) 


z(t)   =  z  =  z  exp(ia)  t 
and       0(t)   =  0  =  0  exp(iu)  t) 

with  the  arbitrary  definitions, 

"z  =  z  exp(-i6) 

°  (7) 

and      0  =  0  exp(-ie) 

Upon  substitution  in  the  original  equations  (1),  a  conscientious  algebraic 
manipulation  yields  the  following  expressions  for  the  "complex"  heaving  and 
pitching  amplitudes, 

-  FS-MQ  -  MP-FR  ,_. 

Z       "       po    rvo  and  9       =       r,o    nD  (8) 


PS-QR  PS-QR 


where, 


(9) 


F     ■     F   (cos  a  -i  sin  a) 
o 

M     =     M   (cos   t  -i  sin  t) 
o 

2 
P     =      (c-aco     )   +  i  b  a) 
v  e  e 

S     =      (C-Aw2   )   +  i  B  w 
e  e 

2 
Q     ■      (g-dw     )   +  i  e  a) 
°  e  e 

R     =      (G-Du)2   )   +  i  E  0) 
e  e 

It  follows  from  (8)  and  (9)  that  the  "complex"  amplitudes  may  be  expressed 
as, 
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/   2   2  Z2 

=  z1  -  i  z_  =  J  z.  +z„   exp[-i  arctan  — ] 


0. 


and  0  =  01  -  i  92  =  7  Q^+Q^     exp[-i  arctan  ~] 
Considering  the  real  parts  of  (6),  we  finally  obtain, 


z  ■  R  z  exp(iu)  t) 
e     c       e 


2  2 

■  Re  J   0n  +0     exp[i(cu  t-arctan  — ] 


ft 


=  0  cos(co  t-e) 
o      e 


(10) 


;2   2  Z2 

z.  +z2   exp[i(co  t-arctan  — ]  (11) 

=  z  cos(co  t-6) 
o      e 

and  0  =  R  0  exp(ia>  t) 


(12) 
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DESCRIPTION  OF  COMPUTER  PROGRAM 
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NOMENCLATURE  OF  COMPUTER  PROGRAMS 


ABAR(I) 

ADDA(I) 

ALPHA 

AVETHE 

BEEB(I) 

BETA 

BMIMAG 

BMNULL 

BMREAL 

BPL 

BSTAR(I) 

CGGC(I) 

CTFST 

CW 

CXFST 

DELTA 

DELV 

DELWL 

DISPL 

DIX(I) 

DMASS(I) 

DRAFT(I) 

DWEIGH(I) 

EMNULL 


see  subroutine  ADMAB 

see  subroutine  COEFF 

phase  lag  of  shearing  force 

average  pitch  amplitude 

see  subroutine  COEFF 

phase  lag  of  bending  moment 

sine  term  of  bending  moment 

bending  moment  amplitude 

cosine  term  of  bending  moment 

ship  (model)  length  between  perpendiculars 

full  beam  (breadth)  of  ship  (model)  station  at  DWL 

see  subroutine  COEFF 

see  subroutine  BENDSH 

wave  celerity 

cosine  term  of  exciting  force 

phase  lag  of  heaving  motion 

increment  of  ship  or  model  speed 

increment  of  wavelength 

ship  (model)  displacement 

d(Sk2k4)/d£ 

station  mass 

station  depth  (draft) 

station  weight 

pitching  moment  amplitude 
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ENOX(I) 

EPSIL 

FNULL 

FROUDE 

GAMMA 

GRAV 

U 

I 

INCRES 

KRIT 

M 

MAXKRI 

MINKRI 

N 

OMEGA 

OMEGAE 

PI 

QUANT (I) 

RO 

SECOE(I) 

SHIMAG 

SHNULL 

SHREAL 

SIGMA 

SKLAM 


station  damping  coefficient 

phase  lag  of  pitching  motion 

heaving  force  amplitude 

Froude  number 

specific  weight  of  water 

gravitational  acceleration 

interval  of  Simpson's  integration 

program  variable 

see  subroutine  BENDSH 

see  subroutine  BENDSH 

no<>  of  ship  or  model  stations 

see  subroutine  BENDSH 

see  subroutine  BENDSH 

highest  station  number 

absolute  wave  frequency 

frequency  of  encounter 

tt  -  3.1415926 


sk2k4 


density  of  water 

section  area  coefficient 

sine  term  of  shearing  force 

shearing  force  amplitude 

cosine  term  of  shearing  force 

phase  lag  of  heaving  force 

see  subroutine  EXCITE  and  BENDSH 
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STFST 

SXFST 

SYMPS 

TAU 

TI(I) 

TIMAG 

TMASS 

TNULL 

TR(I) 

TREAL 

UR(I) 

UI(I) 

V 

VMAX 

VMIN 

WA 

WAVEN 

WL 

XI  (I) 

Y 

YNERT 

ZIMAG 

ZNULL 

ZREAL 


see  subroutine  BENDSH 

sine  term  of  heaving  force 

area  under  any  curve  (see  subroutine  SIMPS) 

phase  lag  of  pitching  moment 

see  subroutine  MOTION 

sine  term  of  pitching  motion 

total  mass  of  ship  or  model 

amplitude  of  pitching  motion 

see  subroutine  MOTION 

cosine  term  of  pitching  motion 

see  subroutine  MOTION 

see  subroutine  MOTION 

ship  or  model  speed 

maximum  speed  of  ship  or  model 

minimum  speed  of  ship  or  model 

wave  amplitude 

wave  number 

wavelength 

longitudinal  coordinate  of  moving  (ship)  system 

see  subroutine  SIMPS 

weight  moment  of  inertia  of  ship  or  model 

sine  term  of  heaving  motion 

amplitude  of  heaving  motion 

cosine  term  of  heaving  motion 
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Description 

(19) 
This  computer  program  was  originally  written  by  K.  Haslum    and  subse- 
quently debugged  and  completed  with  the  aid  of  L.  Vassilopons   .   The  program 
has  the  ability  to  compute  shearing  forces  and  bending  moments  at  any  station 
of  a  ship  or  model  encountering  regular  waves.   Subroutine  BENDSH  which 
essentially  performs  the  above  task,  was  always  bypassed  in  the  main  computa- 
tion cycle  through  the  inclusion  of  a  proper  control  card. 

This  program  is  a  computer  alogorithm  of  the  step-by-step  procedure 
outlined  in  (8)  for  the  calculation  of  ship  (model)  motions  and  bending  moments 
on  the  basis  of  the  Korvin-Krovkovsky  linear  theory  of  pitching  and  heaving. 
Since  the  program  is  extensively  described  in  (19),  only  a  brief  description 
of  its  MAIN  program  and  Subroutines  is  given. 

MAIN  Program 

The  purpose  of  the  MAIN  program  is  to  perform  certain  minor  calculations, 
control  the  computation  cycle  from  subroutine  to  subroutine  and  process  all 
input  and  output  data. 

The  program  output  capabilities  come  under  three  options,  of  which  the 
present  investigation  has  made  sole  use  of  option  No.  1.   The  program  computes 
motions  only  and  assumes  knowledge  of  radius  of  gyration. 

The  main  variable  for  the  program  is  the  wavelength  and  the  basic  computa- 
tion cycle  repeats  the  fundamental  loop,  first  for  different  wavelengths  and 
then  returns  if  the  ship  speed  is  changed.   For  this  investigation  the  ship 
speed  was  zero.   For  the  complete  calculation  the  program  prints  out  the  wave- 
length (WL) 9  ship  speed  (V),  and  frequency  of  encounter  (OMAGAE).   The  units 
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are  consistent  and  correspond  to  those  chosen  in  the  first  READ  statement 
(See  Description  of  Input  Data) . 


Subroutine  ADMAB 

This  subroutine,  obtained  from  the  staff  of  the  Davidson  Laboratory  of 

f  16} 
the  Stevens  Institute  of  Technology,  calculates  according  to  Grim's  theory    , 

a)  the  added  mass  per  unit  length  at  each  station,  which  is  designated 
as  QUANT(I),  and, 

b)  the  ratio  of  the  amplitude  of  the  emitted  wave  to  the  amplitude  of 
heave  at  each  station,  which  is  designated  as  ABAR(I) . 

Although  individually  checked  for  normal  Lewis  ship  sections,  the  range 
of  the  numerical  applicability  of  the  subroutine  has  not  been  thoroughly 
examined. 
Subroutine  COEFF 

This  subroutine  evaluates  numerically  the  expressions  (5)  for  the  coeffi- 
cients of  the  equations  of  motion „   The  correspondence  between  coefficients 
and  program  variables  is  shown  below: 


a  b  c 
d  e  g 
D  E  G 
ABC 


ADDA(l) 
ADDA (2) 
ADDA(3) 
ADDA (4) 


BEEB(l) 
BEEB(2) 
BEEB(3) 
BEEB(4) 


CGGC(l) 
CGGC(2) 
CGGC(3) 
CGGC(4) 


The  damping  coefficient  at  each  station,  denoted  by  ENOX(I)  is  obtained 


from: 


N(£) 


2  —  2 

p£mz 


to 


ee 
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where  A  is  the  quantity  designated  in  subroutine  ADMAG  as  ABAR(I)  and  the 
remaining  symbols  have  their  usual  meaning.  The  numerical  integrations  are 
performed  by  subroutine  SIMPS  with  an  integration  interval  H  =  DX1. 
Subroutine  EXCITE 

Subroutine  EXCITE  computes  the  excitation  force  and  moment  as  discussed 
in  Appendix  B.   The  elemental  excitation  force  at  a  given  ship  station  is 
computed  from  equation  (3)  by  splitting  the  above  expression  into  two  parts. 
The  first  designated  CXFST(I)  is  in  phase  with  cos  co  t  and  is  given  by  the 
first  bracketed  term  of  equation  (3).   The  second  designates  SXFST(I)  is  in 
phase  with  sin  oo  t  and  is  given  by  the  second  bracketed  term  of  equation  (3). 
The  quantities  <$>.   and  <J>_  of  equation  (3)  are  correspondingly  referred  to  as 
FKLAM  and  SKLAM(I).   SKLAM(I)  is  conveniently  subscripted  for  subsequent 
utilization  in  subroutine  BENDSH  if  the  computation  of  bending  moments  is  also 
required.   The  total  excitation  force  and  moment  are  finally  obtained  from 
equation  (4)  by  integrating  over  the  ship  length  the  appropriate  elemental 
values  using  subroutine  SIMPS  with  H  =  DXI. 
Subroutine  MOTION 

Whenever  the  main  program  calls  this  subroutine,  the  excitation  terms  and 

the  coefficients  of  the  equations  of  motion  (1)  are  already  available  and  hence 

the  heaving  and  pitching  motion  amplitudes  and  phase  angles  can  readily  be 

determined.   Subroutine  MOTION  performs  this  task  by  first  expressing  the  set 

(9)  in  the  following  form: 

F  =  F  +  iF^ 
r     i 

M  -  M  +  iM, 
r     i 

P  -  P  +  iP. 
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Q  =  Qr  +  iQ± 

R  =  R  +  iR, 
r     i 

S  -  S  +  iS. 
r     1 

The  various  products  in  the  set  (8)  are  then  evaluated  as  shown  below: 

MQ     =      (M  Q  -M^Q.)   +  i(M  Q>M.Q   ) 

FS     -      (F  S  -F.S.)   +  i(F  S.+F.S   ) 
r  r     i  1  r  x     x  r 

QR     =      (Q  R  -Q,R.)   +  i(Q  R.+Q.R  ) 
r  r     i  i  ^r  x     x  r 

PS     =      (P  S  -P^SJ  +  i(P  S.+P^  ) 
rr     i  x  r  x     ir 

FR     -      (F  R  -F^RJ   +  i(F  R.+F.R  ) 
rr     ii  r  x     x  r 

MP     =      (MP  -M.P.)   +  i(M  P.+M.P   ) 
rr     x  i  r  x     x  r 

so  that  (8)  can  finally  be  computed  from 

7  -  ZREAL  -  iZIMAG 
0  -  TREAL  -  iTIMAG 

where  ZREAL  =  kC  +   B° 

C  +  D 

ZIMAG  =  Af^f 
C  +  D 

TREAL  m    bc^lm 

C  +  D 

TIMAG  =   E°  "  Fj; 
C     +   DZ 
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and 


A 

= 

M  Q 
r  r 

— 

MiQi 

-  F  S 
r  r 

+ 

F  S 
i  i 

B 

= 

M  Q, 
r  i 

+ 

MiQr 

-  F  S. 
r  i 

- 

F.S 
i  r 

C 

= 

Q  R 
r  r 

- 

QiRi 

-PS 
r  r 

+ 

P.S. 

1  X 

D 

= 

Q  R. 
r  1 

+ 

Vr 

-PS4 
r  I 

- 

PiSr 

£ 

- 

F  R 
r  r 

- 

FiRi 

-MP 
r  r 

+ 

M1P1 

F 

= 

F.R 
1  r 

+ 

F  R. 
r  1 

-MP. 
r  i 

- 

M.P 
l  r 

The  amplitudes  and  phase  angles  of  heaving  and  pitching  motions  are  finally 
computed  with  equations  (11)  and  (12),  i.e., 


(ZREAL)2  +  (ZIMAG)2 


0 


(TREAL)2  +  (TIMAG)2 


and 


r   ZIMAG, 

=  arc  tan  {"  ^£aT} 

r      TIMAG , 

=  arc  tan  {"  ^EaI} 


The  correspondence  between  program  variables  and  the  quantities  referred 
to  above  is  as  follows: 


R. 


M    M 


TR(1) 

TI(1) 

TR(2) 

TI(2) 

TR(3) 

TI(3) 

TR(4) 

TI(4) 

TR(5) 

TI(5) 

TR(6) 

TI(6) 
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and 


E  F 
C  D 
A    B  i 


UR(3)  UI(3) 
UR(4)  UI(4) 
UR(5)   UI(5) 


It  should  be  noted  that  TR(2) ,  TI(2),  TR(5)  and  TI(5)  are  furnished  by  sub- 
routine EXCITE. 
Subroutine  BENDSH 

For  a  complete  description  of  this  subroutine  the  interested  reader  is 
referred  to  (19),  which  presents  the  basic  analytical  treatment  involved  in 
the  computation  of  ship  bending  moments  in  regular  waves.   Since  this  subroutine 
is  at  present  not  used  by  the  computer  program  under  description,  we  shall 
refrain  from  further  details. 
Subroutine  SIMPS 

This  subroutine  is  an  integrator  and  calculates  numerically  various 
integrals  and  areas  under  curves  by  using  Simpson's  rule.  At  any  stage  of  the 
computation  cycle,  the  number  of  ordinates  furnished  to  this  subroutine  must  be 
odd.   Thus  any  integral  is  evaluated  in  accord  with: 


f(x)dx  = 


K 


[Y1  +  4Y2  +  2Y3  + 


+  4Y  +  Y 

N    N+1J 


where  Y(I)  are  the  integrand  ordinates  and  I  rises  from  1  to  J.   Y(I),  H  and 
J  must  be  preset  appropriately  to  ensure  that  J  is  an  odd  number. 
Description  of  Input  Data 

The  MAIN  program  includes  six  READ  statements,  numbered  1001  -  1006, 
whose  purpose  is  to  (a)  transmit  the  basic  input  information  required  for  the 
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computations  and  (b)  incorporate  the  control  cards  which  decide  on  the  path  of 
the  computation  and  subsequently  the  form  in  which  the  output  information  is 
to  be  given. 

Unless  specifically  stated  in  the  following,  numerical  data  must  be  of 
the  floating  point  type.   Any  consistent  system  of  units  may  be  used,  the 
decision  being  made  by  the  user  for  the  input  variables  of  the  first  READ  state- 
ment (No.  1001).   At  the  end  of  the  FORTRAN  or  binary  deck  and  after  the  *DATA 
card,  the  basic  information  for  the  computations  should  be  given  in  the  follow- 
ing order,  one  card  per  READ  statement,  unless  otherwise  stated: 

a)  READ  statement  1001,  (I10,4F10.4) 

1)  N       =  number  of  stations  the  ship  is  divided  into. 

This  number  is  also  equal  to  the  highest 
station  number  when  the  latter  runs  from  0 
at  the  FoPo  through  N  at  the  Ac  P.   N  must  be 
an  integer  and  also,  due  to  the  structure  of 
subroutine  SIMPS,  it  must  be  even.   Unless  the 
DIMENSION  statements  are  increased,  N  should 
always  be  less  than  or  equal  to  20 . 

2)  BPL      =  length  between  perpendiculars  of  ship  or  model. 

3)  GAMMA    =  specific  weight  of  water 

4)  GRAV     ■  gravitational  acceleration 

5)  DISPL    =  displacement  of  ship  or  model 

b)  READ  statement  1002,  (3F10.4) 

This  statement  requires  a  number  of  cards  equal  to  N  +  1  giving,  for 
each  station  starting  from  the  F.P.,  the  following  information  per  card: 
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1)  BSTAR(I)   =  full  beam  (breadth)  at  that  station 

2)  SECOE(I)   =  sectional  area  at  that  station  divided  by  the 

area  of  the  circumscribing  rectangle  at  that 
station c 

3)  DRAFT(I)  =  actual  depth  (draft)  of  section  at  that  station. 

c)  READ  statement  1003,  (2F20.4) 

This  card  must  always  be  included  since  it  is  a  control  card  which  dictates 
the  computer  to  bypass  subroutine  BENDSH.   The  following  information  is  assumed 
known  and  must,  therefore,  be  given: 

1)  RADGYR  ■  radius  of  gyration  of  ship  or  model 

2)  XI (I)   =  distance  of  center  of  gravity  of  ship  or  model  from 

the  F.P. 
To  compute  bending  moment,  use  a  blank  card  for  this  READ  statement. 

d)  READ  statement  1004 

This  statement  is  never  encountered  in  the  computations  described 
herein;  therefore  no  input  card  should  be  included  for  this  statement. 

e)  READ  statement  1005,  (3110) 

The  input  card  to  this  statement  must  always  be  included  since  the 
printed  information  serves  as  a  controller  to  the  program.   Ordinarily,  MINKRI, 
MAXCRI  and  INCRES  are  input  variables  for  subroutine  BENDSH  which,  however,  is 
never  used  in  the  present  computations.   For  this  reason,  the  following 
numerical  values  must  be  given  so  as  to  instruct  the  computer  to  bypass  BENDSH 
and  avoid  related  calculations: 

MINKRI   =   -1 

MAXKRI   =    1 

INCRES   =    1 
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f)  READ  statement  1006,  (7F10.4) 

1)  WA     ■  wave  amplitude, 

2)  SWL    =  shortest  wave  length, 

3)  BWL    ■  longest  wave  length, 

4)  DELWL  ■  increment  in  wave  length  from  SWL  to  BWL, 

5)  VMIN   ■  slowest  speed,  ^ 

I   (units  consistent  with  READ  statement 
r   1001) . 

6)  VMAx   ■  highest  speed,   J 

7)  DELV   -  increment  in  speed  from  VMIN  to  VMAX. 

Description  of  Output  Data 

The  program  prints  out  the  absolute  heave  and  pitch  amplitudes,  the  non- 
dimensional  amplitudes,  the  cosine  term,  the  sine  term,  and  the  phase  angles. 
The  particular  modification  of  this  program  used,  computes  relative  bow  motion 
(REBM) .   Only  the  heaving  and  pitching  non-dimensionalized  amplitudes  were 
utilized  in  this  investigation. 


-65- 


Input  and  Output  Data 


Catamaran  Model   -  Directly  Ahead  Seas 

Input 


EXECUTION 

20 

4.8000 

62.4000    32.2000 

48.8050 

0. 

0. 

0. 

•  5200 

.7625 

.1875 

.6156 

.7940 

.2625 

.6750 

.8290 

.2781 

.7000 

.8350 

.3000 

.7000 

.8350 

.3000 

.7000 

.8350 

.3000 

.7000 

.8350 

.3000 

•  7000 

.8350 

.3000 

.7000 

.8350 

•  3000 

•  7000 

.8350 

.3000 

.7000 

.8350 

.3000 

.7000 

.8350 

.3000 

.7000 

.8350 

.3000 

.7000 

.8350 

.3000 

•  7000 

.8350 

.3000 

•  7000 

.8350 

.3000 

.6750 

.8290 

.2781 

•  6156 

.7940 

.2625 

•  5200 

.7625 

.1875 

0. 

0. 

0. 

ENTER  RADGYR 

AND  OTHER 

DATA  NOW 

.9065 

2.4000 

0 

0 

0 

.1250 

.7500 

10.0000      .2500 

0.         1 

0. 
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Listing  of  Program  and  Subroutines 

THE  DYNAMIC  BENDING  MOMENT  IN  REGULAR  WAVES 

DIMENSION  Y(21) ,TR(6) ,TI(6) ,ADDA(4) ,BEEB(4) ,CGGC(4) ,UR(6) ,UI(6) , 
1DMASS(21) ,QUANT(21) ,SKLAM(21) ,BSTAR(21) ,CXFST(21),SXFST(21) , 
2CTFST(21) ,STFST(21) ,XI(21) ,DIX(21) ,EN0XI(21) ,DRAFT(21) ,DWEIGH(21) , 
3SECOE(21),ABAR(21) 

COMMON  Y,J,SYMPS,DXI,ADDA,BEEB,CGGC,ZREAL,ZIMAG,TREAL,TIMAG,ZNULL, 
1TNULL,DELTA,EPSIL,TR,TI,V,DMASS, QUANT, OMEGAE,SKLAM,KRIT,RO,GRAV, 
2BSTAR,CXFST,SXFST, ALPHA, SHNULL, XI, BETA, BMNULL, GAMMA, DIX,M,WA,WAVEN 
3, CW,ENOXI, SIGMA, TAU , FNULL , EMNULL , DRAFT, DWEIGH,SECOE,TMASS,N,UR,UI, 
4ABAR , PI , SHREAL , SHIMAG , BMREAL , BMIMAG , YNERT , BPL 

900  FORMAT ( 110, 4F10. 4) 

901  FORMAT (3F 10. 4) 

902  FORMAT(F10.4) 
CONTINUE 

904  FORMAT (4F2 0.3) 

905  FORMAT (2F20. 4) 

906  FORMAT(7F10.4) 

907  FORMAT(3110) 

922  FORMAT(120) 

903  FORMAT(F9.4,2F10.4,3F9.4,F8.4) 

909  FORMAT (55H  FROUDE  WL/BPL  OMEGAE  HEAVE  PITCH  BEND) 

923  FORMAT(64xI8) 

SHIP-AND  WATER-CHARACTERISTICS 
PI=3, 1415926 

1001  READ  900, N, BPL, GAMMA, GRAV,DISPL 
PRINT  900,N,BPL,GAMMA,GRAV,DISPL 
RO= GAMMA/ GRAV 

FN=N 

DXI=BPL/FN 

M=N+1 

1002  READ  901,(BSTAR(I),SECOE(I),DRAFT(I),I=I,M) 
PRINT  901,(BSTAR(I),SECOE(I),DRAFT(I),I=1,M) 
J=M 

TMASS=DISPL/GRAV 
916  FORMAT (32H  ENTER  RADGYR  AND  OTHER  DATA  NOW) 
1013  PRINT  916 

1003  READ  905,  RADGYR, XI (1) 
PRINT  905,  RADGYR,XI(1) 
IF (RADGYR)  2,2,5 

2  DO  3  I=1,M 

1004  READ  902,  DWEIGH(I) 
PRINT  902,  DWEIGH(I) 

3  DMASS<I)=DWEIGH(I)/GRAV 
HOMENT=0.0 

DO  4  I=1,M 

L=I-1 

FL=L 

4  HOMENT=HOMENT+DWEIGH (I) *FL 
XI  ( 1) =DXI*HOMENT/DISPL 
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5  DO  6  1=2, M 
L+I-l 
FL=L 

6  XI(I)=XI(1)-DXI*FL 
IF(RADGYR)7,7,9 

7  YNERT=0o0 
DO  8  1=1, M 

3  YNERT=YNERT+DWEXGH(I)*(XI(I)*XI(I)+(DXI*DXI(12.)) 
9  YNERT=DISPL*RADGYR*RADGYR 

1005  READ  907,  MINKRI,MAXKRI,INCRES 
PRINT  907,  MINKRI,MAXKRI,INCRES 

1006  READ  906,  WA,SWL,BWL,DELWL,VMIN,VMAX,DELV 
PRINT  906,  WA,SWL,BWL,DELWL,VMIN,VMAX,DELV 
PRINT  909 

V=VMIN 

11  WL=SWL 

12  IF(MINKRI)  14,14,13 

13  KRIT=MINKRI 

14  WAVEN=2.*PI/WL 
CW=SQRTF (GRAV /WAVEN) 
OMEGAE=WAVEN*ABSF  (CVH-V) 
CALL  ADMAB 

CALL  COEFF 

CALL  EXCITE 

FTWQ=-TI(2) 

EMTWO=-TI(5) 

CALL  MOTION 

L=(N/2)+l 

DELTA=DELTA*57 .295779 

EPSIL=EPSIL*57 . 295779 

TNULL=TNULL/ (WA*WAVEN) 

ZNULL=ZNULL/WA 

IF(MINKRI)  16,16,15 

15  CALL  BENDSH 
921  FORMAT ( 7 2H 

1        STATION) 

PRINT  921 

PRINT  923,  KRIT 

SHNULL=SHNULL*BPL/ (GAMMA* (BSTAR(L) *WA) **2) 

BMNULL=BMNULL/ (GAMMA*WA*BSTAR(L) *BPL*BPL) 

ALPHA=ALPHA*57 . 295779 

BETA=BETA*57 . 295779 

16  DWL=WL/BPL 
FROUDE=V/SQRTF (GRAV*BPL) 
CONTINUE 

PRINT  903,  FROUDE,DWL,OMEGAE,ZNULL,TNULL,BMNULL 
IF(MINKRI)  19,19,17 

17  IF(KRIT-MAXKRI)  18,19,19 

18  KRIT=KRIT+INCRES 
GO  TO  15 
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19  IF(WL-BWL)  20,21,21 

20  WL=WL+DELWL 
TO  TO  12 

21  IF(ABSF(V)-ABSF(VMAX))  22,25,25 

22  V=V+DELV 
GO  TO  11 

25  CONTINUE 
GO  TO  1013 
END(1, 0,0, 0,0, 0,0, 0,0, 1,0, 0,0, 0,0,) 
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SUBROUTINE    COEFF  04/18 

SUBROUTINE    COEFF 

DIMENSION    Y(21),TR(6),T1(6)  ,  ADDA(4)  ,  B    ElU  4  )  ,  C3:,  .  (  4  )  ,  UK(      )  ,  U  I  (  6  >  , 
1DMASS(21 ), QUANT ( 21 1 30  )  , SKL4M ( 2 1 ) ,  BSi  A     ( 2  I ) , CXFSI ( 21 ) , SXI  ST(     1)  , 
2CTFST(21),STFST(21),XI(21),DIX(.  II,    N  >XI  (21) lDRAFT(2Ilf DWEIGHI21 > , 
3SECOEU1  ),  ABARI21,  30)  ,WL(3    ),V(.0),hR    QL  (  21  ,  30)  ,  ABARNU     »30), 
4ADMASN(21,30) 

COMMON    Y,  J,SYMPS,DXI  , ADDA,     cEB,CGGC,  Z    E  M  ,  Z 1  MAT,,  TKE  M  ,  T  1  MAG,  ZNULL, 
1TNULL,DELTA,EPSIL,TR,TI, V,  I  MASS,  QUA  •<  T  ,  OMEGAE  ,  SKLAM,  KR  I  F ,RO,  ,RAV, 
2BSTAR,CXFST,SXFST,ALPHA,SH\IULL,XI,iJL  r ., ,  BMNULL  ,  G  \MMA  ,  D  IX  ,  M,  WA  ,  WAVEN 
3,CW,EN0XI  , SIGMA, TAU,FNULL,-MNULL,DK  »F     ,  OWE  IGH,  Sl:  COt ,  F  MA  S  S  ,  N,  UR  ,  J I  , 
4ABAR,PI,SHREAL,SHIMAG,BMREAL,BMIMAG,  Y  IERT  ,  BPL  ,  F~tEQL  ,  AB4KN  , 
5A0MASN,NSTA,K 

SMALL    A    =    ADDAd  ) 

J  =  M 

00    10    1=1, M 
10    Yd)  =  QUANT(I,K) 

CALL    SIMPS 

ADMAS=SYMPS 

ADDA(1)=ADMAS*TMASS 

CAPITAL    A    =    ADDA(4) 

DO    21     1=1, M 
21    Yd  )  =  QUANTd,K)»(XI(  I  )»»2) 

CALL    SIMPS 

ADDA(4)=SYMPS+YNERT/GRAV 

SMALL    0    =    ADDA13)     =    CAPITA!      D    =    ADD;!     ) 

DO    30    1=1, M 
30    Y(I)=QUANT(I,K)«XI(I) 

CALL    SIMPS 

ADDA(2)=SYMPS 

ADDA(3)=ADDA(2) 

SMALL    B    =    BEEB(l) 

DO    40    1=1, M 

ENOXK I)  =  (GAMMA«GRAV»(ABAR(  I,K)»»2))/iUMEGAE»»j) 
40    Yd)  =  ENOXI( I) 

CALL    SIMPS 

BEEB(1)=SYMPS 

CAPITAL  B  =  BEEBI4) 

DO  50  1=1, M 
50  Yd  )  =  ENOXId  )*(XI(  I)»»2) 

CALL  SIMPS 

TEMPORARILY 

BEEB(4)=SYMPS 

DIX(1)=-QUANT(2,K)/(2.«DXI) 

DIX(M)=QUANT(N,K)/(2.»DXI  ) 

DO  55  1  =  2, N 

55  DIXd )  =  ( QUANT ( I- 1 ,K ) -QUANT {  I ♦ 1 , < ) ) / ( 2  -  « i J X  I  ) 
DO  56  1  =  1, M 

56  Yd  )  =  DIX(  I  )  »XI  (  I  )»X1  {  I  ) 
CALL  SIMPS 

BEEB(4)=BEEB(4)-ABSF(V)*(2.»ADDA( 2)  ♦ SYMPS) 
CAPITAL  E  =  BEEB(3) 

DO  60  1  =  1, M 
60  Y{ I )=DIX{  I  )*XI(  I  ) 
CALL  SIMPS 

ETHRE=SYMPS»(ABSF(V)  ) 
DO  63  1=1, M 
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63  Yd  >=ENOXI(  I  )»XI  (  I  ) 
CALL  SIMPS 
E0NE=SYMPS 

ETW0  =  2.MA6SF(V)  )»ADMAS 
BEEB(3)=E0NE-ETHRE 
SMALL  E  =  BEEB(2> 
BEEB(2l=BEEB(3)-ETWO 
SMALL  C  =CGGC(I) 
DO  70  1=1, M 

70  Y(I»=BSTAR(I) 
CALL  SIMPS 
CGGC(1)=GAMMA»SYMPS 
CAPITAL  C  =  CGGC(A) 
DO  71  I=1,M 

71  Yd  )=BSTAR(  I  )»(XI(  I  )««2) 
CALL  SIMPS 

CGGCU>=GAMMA«SYMPS-<  ABSFU)  »»HQEB(  ) 
CAPITAL  G  =  CGGC(3) 

DO  80  1  =  1, M 
80  Yd  »=BSTAR(I)*XI(  I) 
CALL  SIMPS 
CGGC(3)=SYMPS»GAMMA 
SMALL  G  =  CGGC(2) 

CGGC(2)=CGGC(3»-(ABSF(V» )*BEEB(1I 
RETURN 
END (1,1,0, 0,0,0, 1,1,0,1,0,  ,0,0,0) 


JOB  TIME  =    0.37  MIN. 
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SUBROUTINE    EXCITE  Q<*/1 

SUBROUTINE    EXCITE 

DIMENSION    Y(21),TR(6),TI(6) , ADDA(4) , B    EB 14) tCGGC (4) , UR(     )  ,UM6), 
1DMASS(21  )  ,QUANT(2i,30),SKL  vM(21) ,BSI A     (  2  1  )  ,  C  XFS  I  (  2  i  J  ,  SXF  ST  (  Z.  1  )  , 
2CTFST(21  )  ,STFST(21)  ,XI(  21),DIXdl),    NTXK21)  ,  DRAFTl  21 )  ,  DHEIGH(  21 )  i 
3SECOEI21  >,ABAR(21,30)  , WL { 3  J  )  , V ( 20  )  , F  R    QL ( 2 1 , 30 ) , A6ARN < 2     ,30)  , 
4ADMASN(21,30) 

COMMON    Y, J,SYMPS,DXI , ADDA, -eEB,CGGC , Z    EAL.ZIMA:,, TRE  AL  ,  1"  I  MAG,  ZNULL  , 
lTNULL,DELTA,EPSIL,TR,TI,V,t>MASS,QUA  .  T , OMEGAS t SKI  AM, KR I  X , RO, GRAV , 
2BSTAR,CXFST,SXFST,ALPHA,SHWJLL,XI,iJ    T  v,  BMNULL  ,  G  \MMA,  D  I  X  ,  M,  W  \ ,  WAVEN 
3,CW,EN0XI,SIGMA,TAU,FNULL,=MNULL,0K\Fi , DWEIGH, SECOE, TMA SS.N, UR , UI , 
4ABAR,PI,SHREAL,SHIMAG,BMREAL,BMIMAG, Y  iEKT,BPL,F    CQL,ABA?>N, 
5ADMASN,NSTA,K 

J  =  M 

DO    90    1=1, M 

FKLAM=( GAMMA »BSTAR( I  )-( (WAVEN»CW)»«     )  » QUANT (  I  ,  K )  ) *WA 

SKLAM( I)=ENOXI( I ) -( D I X{ I ) »\BSF ( V ) ) 

CXFSTt  I  )=(FKLAM*SINF( WAVEN*XI(  I)  )  ♦(  W  At/EN«CW*WA  )  «SKLAM(  I  ) 
1*C0SF(WAVEN«XI( I ) ) ) • ( EXPF ( -W AVEN*DR AF i  ( I  I »SECOc  {  I ) )  ) 

SXFST( I)=(FKLAM»COSF(WAVEN«XI( I )  )-( WAVEN«CW«WA ) « SKLAM  I ) 
1«SINF(WAVEN»XI( I ) ) )* (EXP Ft- WAVEN* DR^F  r ( I ) «SEC06 ( I ) ) ) 
90    Yd  )=CXFSTd  ) 

CALL  SIMPS 

FONE 

TR(2)=SYMPS 

DO  92  1=1, M 

92  Yd  >=SXFST(  I  ) 
CALL  SIMPS 
-FTWO 

TI(2)=-SYMPS 
DO  93  1=1, M 

93  Y( I )=CXFST(  I  )  »X1  (  I  ) 
CALL  SIMPS 

EMONE 

TR(5)=SYMPS 
DO  94  1=1, M 

94  Yd  )=SXFST(  I  )»XI  (  I) 
CALL  SIMPS 

-EMTWO 

TI (5)=-SYMPS 

SIGMA=ATANF(-TI(2)/TR(2)  ) 

TAU  =  ATANF(-TI(5)/TR(5)  ) 

FNULL=SQRTF(TR(2)«»2+TI{2)*«2> 

EMNULL=SQRTF(TI(5)**2+TR(5)»»2) 

RETURN 

END( 1,1, 0,0,0,0, 1,1,0,1,0,  ,0,0,0) 


JOB  TIME  =   01.20  MIN. 
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SUBROUTINE  MOTION  0^/13 

SUBROUTINE  MOTION 

DIMENSION  Y(21),TR(6) ,TI (6) , ADUA(4)  ,B  Eb ( 4) , CGG„ ( 4 ) » UR(  ),UI(6), 
1DMASS(21 ) ,QUANT( 21,30) fSKL<M( 21 )  ,BS I  A • U 1 ) , CXFST ( 21 ) , SXF ST(  < L  I  , 
2CTFST(21)iSTFST(21)iXI(211f  DIX(Zllt~NuXI  (21ltDi4AFTC;  I,  f>W£IGH(2i  )  t 
3SECOE(21 ) ,ABAR(21,30),WL(3  >,V120),FR  QL ( 21 , 30) , ABAKNC2  ,30), 
4ADMASN(21,30) 

COMMON  Y, J,SYMPS,DXI , ADDA,  <CEB, CGGL ,  Z^ EALt Z J  MAG ,  1  REAL , T I  MAG,  /NULL, 
1TNULL, DELTA,  EPSI L VTR, TI , V, "MASS, QUA ^ T , OMEGA r ,  SKL.AM,  KRIT  ,  RQ, GRAV, 
2BSTARtCXFST,SXFST,ALPHA, SHMULL , X  I , rt _ T  ,  rtMNUl L , G aMMA, DIX , M, W4, HAVEN 
3,CW,EN0XI,SIGMA,TAU,  FNULL  ,  lMNULL  ,  OK  F  i  ,  I) WE  I  GH  ,  1>  -COE  ,  TM*  r,  S  ,  N,  UR  ,  U  I  , 
4ABAR,PI ,SHREAL,SHIMAG,BMREAL,BMIMAG, Y  ,ERT,BPL,F  <EQL,ABARN, 
5ADMASN,NSTA 

DO  105  J  =  l,4 

GO  TO  1100, 101, 102, 103),  J 

100  1=1 
GO  TO  104 

101  1=3 
GO  TO  104 

102  1=4 
GO  TO  104 

103  1=6 

104  TR( I)=CGGC( J) -ADDA ( J)*(OME GAE*»2  ) 

105  TIU)=BEEB(J)«OMEGAE 
DO  110  1=2,3 
DO  110  K=l,2 
IF(K»I-4>  108,110,108 

108  IPK=I+K 

UR(  IPK)=TR(I)»TR(K  +  3)-TI(  I  )»TI(K*3)  ■  T  (K)»TR(  I^)*l  I  (K)  *  I  K  1*3) 
Ul(IPK)=TR(I)«TI(K+3)+TI(I)*TR(K+3)-T(K)»TI(H  )-TI(K)«TK(I+3) 

110  CONTINUE 
DO  111  1=1,6 
UR(I )=UR(I»/10000000000. 

111  UI(I)=UI(I)/10000000000. 
DENUM=UR(4)»»2+UI(4)«»2 
ZREAL=(UR(5)»UR(4)+UI(5)»UI (4) )/DEN>lM 
ZIMAG=(UR(5)»UI(4)-UI(5)»UK(4))/DENUM 
TREAL=(UR(3)»UR(4)+UI(3)»UI (4) )/DENUM 
TIMAG=(UR(3)«UI(4)-UK3)*Uk(4) J/DENUM 
ZNULL=SQRTF(ZREAL**2+ZIMAG»»2) 
TNULL=SQRTF(TREAL**2+TIMAG»«2) 
DELTA=ATANF(ZIMAG/ZREAL) 
EPSIL=ATANF(TIMAG/TREAL) 
RETURN 
END (1,1,0,0,0,0,1,1,0,1,0,0,0,0,0) 

JOB  TIME  =    0.83  MIN. 
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SUBROUTINE  BENDSH  02/02 

SUBROUTINE  BENDSH 

01  HENS  I  ON  Y<21  )  ,  TR  <  6  WT  I  <  6  )  ,  AODA  (  4  )  ,  BEEB  »  4  )  , CGGC < 4 ) , UR ( 6 ) JUI < 6 )  , 
10MASS<21),QUANT(21),SKLAM<2l),BSTAR<2i),CXFST<21),SXFST(21l, 
2CTFST<21) ,STFST(2l ) ,XI(21> ,DIX(2l),ENOXI(21),DRAFT(21),DWPIGH<21) . 
3SECOE<21),ABAR(21),ACREAH21),ACIMAG<21) 

COMMON  Y,  J,SYMPS,DXI,ADDA,BEEB,CGGC,ZREAL,.ZIMAG,TREAL,TIMAG,ZNULL, 
1TNULL,CELTA,EPSIL,TR,TI, V,  DMASS  ,  QUANT  ,  OMEGAE  ,  SKLAM  ,KR  I  T  ,R0  ,-GRAV  , 
2BSTAR,CXF ST, SXFST, ALPHA, SHNULL, XI, BETA, BMNULL , GAMMA , D I  X, MV WA ,WAVEN 
3,CW,ENCXI,SIGMA,TAU,FNULL,EMNULL»DRAFT,DWEIGH,  SECOE, TMASS^N.UR ,UI , 
4ABAR,PI,SHREAL,SHIMAG,BMREAL,BMIMAG,YNERT,BPL 

X*(-1.)»»KRIT 

IF ( XI  111,111,112 

111  J*KRIT 

GO  TO  113 

112  J«KRIT*1 

113  JENS-MRIT+1 

DO    114    1=1, JENS 

ACREAL<n*(ZREAL*TREAL»XI< II )•( 0MEGAE»»2 ) 

ACIMAG(n»(ZIMAG  +  TIMAG*XI(  I)  l»C  0MEGAE»«2  > 

CTFSTU)*ACREAL<  I)»QUANT(I) 
l*(QUANTUI»2.-ABSF<V)«0MEGAE)»TIMAG-SKLAMtI )•<( < Z IMAG+TIMAG* 
2Ximi«0HEGAE)-ABSF(VI»TREAL)-(GAMMA»BSTAR|  III* <ZREAL*TRE*L»XI U  )) 
3+OXFSIU) 

STF$T< I )»ACIMAG<  I  )»QUANT< I  ) 
l-(  QUANT  m«»2.»ABSP(V)*0MEGAE)»TREAL*SKLAMm-<(  (  ZREAL*TREAL» 
2Xl(m«0MEGAE)+A6SF(V)«TIMAG)-(GAMMA*BSTAR(m«(ZIMAG+TIMAG*XI<I  )) 
3*SXFSTU> 

114  Y(I)»GTFSTm 
CALL    SIMPS 
SHREAt*Ok 

DO    200    1=1, J 

200  shrEal»shre al+acreal <  n *dmass< i ) 

SHREAL=SHREAL-(X«ACREAL( J ENS  )* DM ASS ( JENS)/2.  J 
SHREAL*SYMPS-(X»CTFST<KRIT*n»DXl/2.)*SHREAL 

117  00  118  1*1, J 

118  Y(II»StFSTU) 
CALL  SIMPS 
SHlMAG'Ow 

00  250  1=1, J 
250  SHlMAG*SHlMAG*ACIMAG<I)»DMAS$m 

SHIMAG=SHIMAG-X»ACIMAG(JENS)»0MASS(JENS)/2. 
SHIMAG»SYMPS-(X»STFST(KRIT*1I«DXI/2.)»SHIMAG 

121  ALPHA*ATANF(SHIMAG/SHREAL) 
SHNULL*SQATFCSHREAL»«2+SHIMAG»»2) 
00  122  1=1, J 

122  Y(n»CTFST(II«(XI(  I  )-XI  (KRIT*1)  I 
CALL  SIMPS 

BMREAL-O* 

00    300    1*1, J 
300    BMREAL=BMREAL+ACR6AL< I )*DMASS<  I  I  •  I  X  I  f  I  )-XI(  JENS)  > 

8MREAL= BMRE AL+ ACRE AL ( JENS )»DMASS< JENS) »DX 1/8. 

BMREAL=SYMPS*(CTFST<KRIT*1)»0XI»DXI/8.I*BMREAL 
125    00    126    1=1, J 
12  6    Yd  )*STFST(  I )»( X  H  I  )-X H KR IT+1 )  ) 

CALL    SIMPS 

BM[MAG=0. 
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SUBROUTINE    BENDSH  02/02    12 

00    350    1*1, J 
350    BMIMA<3=BHMAG+ACIMAG(  I )  »DMASSH  )•  (  X  I  (  I  >-XI(  JENS)  ) 

BMIMAG  =  BNIMAG-i-ACIMAG(JENS)»DMASS<JENS)»DXI/8. 

BMIMAG»SYMPS*<STFST(KRITU)»DXI»DXI/8.)+BMlMAG 
129    BETA*ATANF(BMIMAG/BMREAL) 

BMNULL=SQRTF<BMREAL»»2+BMIMAG»»2) 

RETURN 

END (  1,0,0, 0,0,0,0, 0,0, 1,0,0,0 ,0,0) 

JOB  TIME  =     .05  MIN. 


SUBROUTINE  SIMPS  0^/18 

SUBROUTINE  SIMPS 

DIMENSION  Y(21),TR(6),TI(6) , ADDA(4) ,H.EB(4I ,CGGC(4) ,UK(  ) ,UI (6) , 
1 DMAS S (  21 ) , QU ANT  (  2 1 , 30 ) , SKL  AM ( 2 1 )  ,  BS  F  A  .  (  2 1 ) , u  XF  S I ( 2  i  > , SX E  S  T ( ^ 1 )  , 
2CTFST(21»,STFST(21) ,XI( 21 I 1 DIXI 21 ) , rNOXI ( 21 » f DkAFTI 1 1 1 , DWEIGHl 21 ) i 
3SEC0E(21),ABAR(21i30) ,WL(3C ),V(20),FRQL(21,30) , ABARNU- ,301 , 
4ADMASN(21,30) 

COMMON  Y, J,SYMPS,DXI , ADDA, BEEB,CGGC , Z ; EAL , 1 1  MAG , TREAL.T IMAG,2NULL, 
1TNULL t DELT At EPSIL i TRtTIf VfDMASSi QUANT, OMEGA tiSKL AM, KK I  X , RO,  GRAV , 
2BSTAR,CXFST,SXFST,  ALPHA,  SHMULLfXI  fBET  < ,  rtMNUI  L  ,  G  \MMA  ,  L)  I  X  ,  M  ,  W  A  ,  WAVEN 
3,CW,EN0XI,SIGMA,TAU,FNULL,lMNULL,DR \Fi  ,DWEIGHfSLCOE, 1  MA S 5 , Nt UR , U 1 , 
4ABAR,PI,SHREAL,SHIMAG,BMREAL,BMIMAG, Y  JER T , BPL , F^EQL , AflAn N, 
5ADMASN,NSTA,K 

SUM=Y(1)+Y(J) 

KARL=J-1 

DO  151  1=2, KARL 

X=(-i.)»»I 

IF(X)  151,150,150 

150  Yd  1=2. •Yd) 

151  SUM=SUM+2.»Y( I) 
SYMPS=(DXI/3. )»SUM 
RETURN 

END( 1,1, 0,0, 0,0, 1,1, 0,1,0,  •, 0,0,0) 

JOB  TIME  =   01.56  MIN. 


-93- 


APPENDIX  D 
MODEL  RADIUS  OF  GYRATION 

The  procedure  used  to  determine  the  radius  of  gyration  of  the  models  is 
described.   No  information  was  available  on  the  full  scale  radius  of  gyration 
for  the  catamaran  or  ALVIN.   Therefore  the  radius  of  gyration  was  assumed  to 
be  0.25  of  the  length.   Because  of  the  small  weight  of  ALVIN  compared  to  weight 
of  instrumentation,  which  had  to  be  positioned  at  the  center  of  gravity,  the 
largest  radius  of  gyration  that  could  be  obtained  was  0.232  of  the  length. 
The  catamaran  model  was  constructed  with  built-in  lead  ballast,  and  even  with 
a  drastic  alteration  of  the  model,  the  largest  radius  of  gyration  was  0.189  of 
the  length. 


2L 


(V>  x  ^ 


(  W  ,;k  ^ 


(W) 
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Equipment; 

2  springs  of  same  stiffness  k 

1  stopwatch 

string,  tape,  set  of  ballast  weights 
Procedure: 

a)  Suspend  the  two  springs  from  the  clamps  which  are  attached  to  convenient 
points  on  the  ceiling.   The  distance  between  the  two  springs  should  be 
2L. 

b)  Make  sure  the  ship  model,  without  ballast,  rests  on  level  keel  (or  at 
desired  trim)  when  put  in  water,  and  determine  the  center  of  gravity 
of  the  modelo   Assume  that  the  center  of  gravity  coincides  with  the 
center  of  flotation.   Now  put  all  ballast  in  model  symmetrically  about 
the  C.G. 

c)  Measure  off  the  distance  L  on  each  side  of  the  center  of  gravity  of  the 
ship  model.   L  can  be  of  any  convenient  length. 

d)  Suspend  the  model  from  the  springs,  using  string  wrapped  around  the 
model  and  secured  in  place  by  tape  as  shown  in  the  figure  above. 

e)  Force  the  model  to  an  up  and  down  motion  (heave)  and  count  with  the 

stopwatch  the  number  of  oscillations  per  minute.   Thus,  determine  f  . 

n 

f)  Force  the  model  to  a  pitching  motion  and  count  with  the  stopwatch  the 

number  of  oscillations  per  minute.   Thus  determine  f  .   It  is  not 

P 

necessary  to  achieve  a  pure  pitching  motion,  completely  uncoupled  from 
heave,  because  the  period  of  heave  does  not  affect  the  period  of  pitch. 

g)  Let  W  be  the  weight  of  the  ship  and  w  the  weight  of  each  of  the  springs. 
The  moment  of  inertia  of  the  system  can  be  approximated  by: 
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I   =   (W+w)k  2   and   I      Wk  2+wL2  (13) 

s  s  S        g 

where  k  and  k  are  the  radius  of  gyration  of  the  system  and  the  model, 
s      g 

respectively . 

2     W+w     2  w   2 
Therefore:        k  Z  =  ™   k  -  g  IT  (14) 

g       W     s   W 

The  radius  of  gyration  of  the  system  is  given  by: 

ks  =  /  L,  since  fR  =  \     &  f  =  \    ^       (15) 
P  s 

where  M  is  the  mass  of  the  system  and  k  the  stiffness  of  each  spring. 

Hence, 

f  2 

kg2  =  [T  0  "  w]l2  (16) 

For  a  given  weight  of  model,  the  heaving  period  f   remains  constant. 

f  H 

H 
So,  we  can  change  the  ~  ratio  by  changing  the  pitching  period.   To 

P 
increase  f  ,  dispose  ballast  weights  closer  to  C.G.,  making  sure  they 

are  still  symmetrically  disposed  about  the  C.G.   To  decrease  f  ,  dispose 

H 
ballast  weights  further  away  from  the  C.G.  When  ~p~     equals  the  value 

P 
given  from  the  above  relation,  the  ballasting  procedure  of  the  ship  is 

over. 
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APPENDIX  E 
SUMMARY  OF  EXPERIMENTAL  DATA 

The  following  data  is  a  summary  of  the  experimental  data  obtained  at  the 
M.I.T.  Ship  Model  Towing  Tank.  It  represents  data  read  from  Sanborn  Oscillo- 
graph tapes. 

The  catamaran  single  hull  and  catamaran  data  are  for  both  directly  ahead 
and  astern  seas.   (Tables  V-VIII) .   Table  IX  is  for  ALVIN  in  directly  ahead 
seas  and  Table  X  for  directly  astern  seas.   Table  XI  is  ALVIN  data  in  directly 
ahead  seas,  while  in  the  recovery  position.   Table  XII  is  ALVIN  pitching  data 
and  the  catamaran  heaving  data,  while  in  the  recovery  position. 
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TABLE  V 


Catamaran  Single  Hull  Experimental  Data,  h  • 

=  0.375  inches 

0  rad) 

X(ft.) 

zo(in) 

0o(deg) 

z  /h 
o  o 

o 

Kh 
o 

3.0 

0.0 

0.125 

0.0 

0.034 

4.0 

0.04 

0.25 

0.1068 

0.089 

4.5 

0.0025 

0.5 

0.0066 

0.200 

6.0 

0.075 

0.5625 

0.2000 

0.300 

10.0 

0.1875 

0.75 

0.5000 

0.705 

TABLE  VI 

Catan 

laran  Experimental  Data, 

h  =  0.375  inches 
o 

A(ft.) 

zo(in) 

eo(deg) 

z  /h 
o  o 

0  (rad) 
o 

Kh 
o 

3.0 

0.0400 

0.25 

0.1068 

0.067 

4.0 

0.0250 

0.80 

0.0667 

0.284 

4.5 

0.0375 

0.75 

0.1000 

0.301 

6.0 

0.0875 

0.625 

0.2333 

0.333 

10.0 

0.3000 

0.825 

0.8000 

0.775 
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TABLE  VII 


Catamaran  Single  Hull 

Experimental 

.  Data,  h 
o 

=  0.75  inches 

9  (rad) 

A(ft.) 

zo(in) 

0o(deg) 

z  /h 
o  o 

o 

Kh 
o 

3.0 

0.0 

0.20 

0.00 

0.027 

4.0 

0.075 

0.75 

0.10 

0.134 

4.5 

0.025 

1.00 

0.334 

0.201 

6.0 

0.1625 

1.69 

0.2165 

0.452 

10.0 

0.3750 

1.375 

0.500 

0.613 

TABLE  VII] 

[ 

Catanu 

iran  Experi 

mental  Data, 

h  =  0.75 
o 

inches 

0  (rad) 

X(ft.) 

zo(in) 

0o(deg) 

z  /h 
o  o 

o 

Kh 
o 

3.0 

0.10 

0.25 

0.1333 

0.0334 

4.0 

0.05 

1.00 

0.0667 

0.1780 

4.5 

0.0625 

1.25 

0.0834 

0.2500 

6.0 

0.200 

2.25 

0.2670 

0.6000 

10.0 

0.600 

2.00 

0.8000 

0.8910 

-99- 


TABLE  IX 


ALVIN  Experimental  Data,  Directly  Ahead  Seas,  h  =  0.375  inches 


X(ft.) 


zQ(in) 


0o(deg) 


z  /h 
o     o 


0   (rad) 
o 

Kh 


1.0 
2.0 
3.0 
5.0 
10.0 


075 

,150 
,150 
,08 


0.5 
0.6 
1.5 
2.7 
5.4 


0.2 

0.4 
0.4 
0.2135 


0.0045 
O0IO68 
0.4000 
1.2000 
4.8000 


IABLE_J£ 


ALVIN  Experimental  Data,  Directly  Astern  Seas,  h  =  0.375  inches 


X(ft.) 


zo(in) 


eQ(deg) 


z  /h 
o  o 


0  (rad) 
o 

Kh 


1.0 
2.0 
3.0 
5.0 
10.0 


0.01 
0.05 
0.06 
0.13 


0.9 
0.7 
1.4 
4.8 


0.0266 
0.133 
0.160 
0.347 


0.08 
0.1245 
0.375 
2.12 


-  indicates  data  not  obtained 
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TABLE   XI 


Recovery  Experimental  Data  Directly  Ahead  Seas,  h  =  0.375  inches 

A(ft) 

ALVIN 
zo(in) 

ALVIN 
0o(deg) 

z  /h 
o  o 

0  (rad) 
o 

Kh 
o 

1.0 

2.0 
3.0 
A. 5 
5.0 
10.0 

0.030 
0.055 
0.100 
0.0A0 
0.005 
0.075 

0.60 
1.00 
A.A0 
1.50 
2.00 
5.00 

0.008 

0.01A7 

0.0267 

0.0107 

0.0133 

0.0200 

0.005 
0.178 
1.175 
0.600 
0.890 
A.A5 

TABLE  XII                             i 

'1 

Recovery  Experimental  Data  Directly  Ahead  Seas,  h  =  0.375  inches 

X(ft.) 

Cat. 
zQ(in) 

ALVIN 
0o(deg.) 

z  /h 
o  o 

ALVIN 

0  (rad) 
o 

Kh 
o 

3.0 
A.O 

A. 5 

6.0 

10.0 

0.050 
0,010 
0.025 
0.125 
0.225 

1.7 
1.0 
0.8 
1.6 
2. A 

0.133 
0.027 
0.067 
0.33A 
0.600 

0.A5A 
0.356 
0.320 
0.853 
2.1A0 
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